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Abstract. We give the path model solution for the cluster algebra variables of the A r T-system with 
generic boundary conditions. The solutions are partition functions of (strongly) non-intersecting paths 
on weighted graphs. The graphs are the same as those constructed for the Q-system in our earlier work, 
and depend on the seed or initial data in terms of which the solutions are given. The weights are "time- 
dependent" where "time" is the extra parameter which distinguishes the T-system from the Q-systcm, 
usually identified as the spectral parameter in the context of representation theory. The path model is 
alternatively described on a graph with non-commutative weights, and cluster mutations are interpreted as 
non-commutative continued fraction rearrangements. As a consequence, the solution is a positive Laurent 
polynomial of the seed data. 



1. Introduction 

In this paper we study solutions of the T-system associated to the Lie algebras A r , which we write in 
the following form: 

(1-1) T a j.k+iT a .j t k-i = T a j + i t kT a j^i t k + T a+ i j,kT a -i j.k, 

where j, k G Z, a 6 I r = {1, r}, and with boundary conditions 
(1-2) To,j,k = T r +i,j,k = 1, j, k e Z. 

We consider these equations to be discrete evolution equations for the commutative variables {T a j t k} in 
the direction of the discrete variable k. 

Originally, this relation appeared as the fusion relation for the commuting transfer matrices of the 
generalized Heisenberg model [TJ [17] associated with a simply-laced Lie algebra g, where it is written in 
the form 

(1-3) T a J^ + lT a J,k-l = ^aj + liTttj-li: — Tfij^k ! 

with appropriate boundary conditions. The matrix C is the symmetric Cartan matrix of one of the Lie 
algebra of type ADE. Our relation (jl.ip is obtained by a rescaling of the variables T a j t k and specializing 
to the Cartan matrix of A r . 

With special initial condition at k = 0, it has been proved that the solutions to (|1.3p are the g-characters 
[TO] of the Kirillov-Reshetikhin modules of the affine Lie algebra U q (sl r +i) [T9] , 

The T-system also appears in several other contexts. Of particular relevance here is the fact [18] that 
the system is a discrete integrable equation, the discrete Hirota equation. It is therefore to be expected 
that the system has a complete set of integrals of motion, and that it is exactly solvable. This equation 
also appears in a related combinatorial context, as the octahedron equation, which was studied by [161 121] . 

In this paper, we do not impose any special boundary conditions, but express the general solution of 
the T-system in terms of arbitrary initial conditions. For example, initial conditions can be chosen by 
specifying the values of the parameters T a j ^ at k — and k = 1, or a more exotic boundary can be 
specified. To solve the system, we use a path model which is a simple generalization of the path model we 
constructed for the solutions of the Q- system of A r [5] [3] . 
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In our previous work, we constructed a set of path models, and proved that the solutions of the Q-svstem 
of A r [15], 

Qa,k+lQa,k-l = Qa,k + Qa-l,kQa+l,k, Qo,k = Qr+l,k = ljfc 6 Z, Of G I r , 

are the generating functions for paths on a positively weighted graph, where the weights are a function of 
the initial conditions. 

With special initial conditions at k — and k = 1 (together with a rescaling as in (| 1 . 3|) which restores 
the minus sign in the second term on the right hand side of the Q-system) , the solutions are the characters 
the finite- dimensional, irreducible modules of A r with highest weights which are multiples of one of the 
fundamental weights. 

Note that this Q-svste 111 is obtained by "forgetting" the spectral parameter j in Equation Thus 
the T-system can be regarded as an afhnization or g-deformation of the Q-system, and the path model we 
present here is therefore a deformation of the path model for the Q-system. 

Without fixing any special initial conditions, it was shown in [13j that the solutions of the Q -S y s tem 
are cluster variables in a cluster algebra |8j. We showed in [4] that all Q-systems, corresponding to any 
simple Lie algebra, can be formulated as cluster algebras. Thus, the solution of the Q-system in terms of 
the statistical model allowed us to prove the positivity conjecture of [5] for these cluster variables. In fact, 
as we showed in [6], the solutions are related to the totally positive matrices of |9] corresponding to pairs 
of coxeter elements. 

Similarly, we showed in {4j that a large class of equations which we call generalized bipartite T-systems 
can be formulated as cluster algebras. Equation is perhaps the simplest example of such a system. 
Motivated by our statistical model introduced in [5] , we introduce a path model which provides us with the 
solution to the T-system, in terms of a set of initial conditions, as the partition function of a path model 
with time-dependent (or non-commutative) weights. Here, we refer to the variable normally identified as 
the spectral parameter as the time parameter, as it is a natural interpretation from the point of view of 
paths. 

This paper is organized as follows. In Section^ we review the necessary definition of a cluster algebra. 
We recall our formulation [4] of T-systems as cluster algebras. We describe the conserved quantities of 
the T-systcm in terms of discrete Wronskian determinants in Section [3] We define a generalized notion 
of hard particle models on a graph in Section [4] and identify the conserved quantities as hard particle 
partition functions on a specific graph. In Section [SJ we use our conserved quantities to write the solutions 
of the T-system as the partition functions of paths on a weighted graph. The weight of a step in a path 
depends on the order in which the steps are taken, that is, the weights are time-dependent. The solutions 
are written as functions of the fundamental initial data, and the graph is the same as the one used in the 
Q-system solution. Positivity of the T-system solutions in terms of the fundamental seed variables follows 
from this formulation. 

To prove the positivity in terms of other seeds, we give a formulation of our model in terms of non- 
commutative weights in Section [HI We are then able to describe the solutions of the T-system as a function 
of other seed data as partition functions on new graphs with weights which depend on the mutated seeds. 
The key to the construction is an operator version of the fraction rearrangement lemmas used in [4] . These 
rearrangements are equivalent to mutations in the case of the Q-system. Here, they are equivalent to 
compound mutations. We are thus able to write the T-system solution explicitly in terms of its initial 
data, for a subset of cluster seeds. 

This paper should be considered as a (special case of) non-commutative generalization of our work 
on the solutions of Q-system [5J [6]. In particular, the graphs on which we build our path models are 
the same as for the Q-system, and the only difference is the time-dependence or non-commutativity of 
the weights. The various key properties, such as the rearrangement lemmas for continued fractions and 
the generalization of the Lindstrom-Gessel-Viennot theorem for strongly non-intersecting paths, all have 
straightforward non-commutative counterparts which are used here. 
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2. T-SYSTEMS AS CLUSTER ALGEBRAS 

2.1. Cluster algebras. We use the following definition of a cluster algebra [El[22]> slightly specialized to 
suit our needs in this paper. 

Let S C S be two discrete sets (possibly infinite) and consider the field J of rational functions over Q 
in a set of independent variables indexed by S. 

We define a seed in 5" to be a pair (x, B), where x = {x m : m £ S} is a set of commuting variables, and 
B is an integer matrix, with rows indexed by S and columns indexed by S. The matrix B, which is the 
square submatrix of B made up of the rows of B indexed by S, is skew symmetric. 

The cluster of the seed (x, B) is the set of variables {x m : m £ S}, and the coefficients are the set of 
variables {x m : m £ S \ S} . 

Next, we define a seed mutation. For any m £ S, a mutation in the direction m, \x m : (x,B) <— > (x',B'), 
is a discrete evolution of the seed. Explicitly, 

• The mutation /x m leaves x n with n ^ m invariant, and updates the variable x m only, via the 
exchange relation 

(2.i) .<•:„ .<•„/ (i i •<•;," • n ; ; 

Vugs ties 

where [n] + = max(n,0). 

• The exchange matrix B' has entries 



(2-2) = 



—Bi j if i = m or j = m; 

Bi.j + s\gn(Bi, m )[Bi, m B m ^] + ) otherwise. 



Note that we only define mutations for the set S, and not for the coefficient set S\S. That is, coefficients 
do not evolve. 

Fix a seed (x, B) and consider the orbit 1 C J of the cluster variables under all combinations of the 
mutations /i m , m £ S. The cluster algebra is the Z^ 1 ]- subalgebra of 1 generated by X, where c is the 
common coefficient set of the orbit of the seed. 

Remark 2.1. The particular system which we solve in this paper does not require us to have a coefficient 
set, that is, we can set S = S. However, to make more direct contact with representation theory, it is 
desirable to have the coefficient set be enumerated by the roots of the Lie algebra. In this context, we need 
to set the values of the coefficients to the special points —1. 

Cluster algebras can be considered to be discrete dynamical systems, which is the point of view we adopt 
in this paper. 

2.2. bipartite T-systems as cluster algebras. In this section we review some of the definitions of 
Appendix B of 4J, where generalized bipartite T-systems were shown to have a cluster algebra structure. 

Definition 2.2. A generalized bipartite T-system is a recursion relation for the commuting, invertible 
variables {T a j-k}, where a £ I r and j, k £ Z, of the form 

(2.3) T a j-k+iT a j-k-i — T a j + i^T a j^i-k + q a J^[(7a' °' ,Q 
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where A is an incidence matrix, that is, a symmetric matrix with positive integer entries. 

The matrix A is generally of infinite size, unless special boundary conditions are imposed on the system 
which truncate the range of the variables j. We do not impose such boundary conditions in this paper, 
although they are clearly of interest [12, 20 . The symmetry of A is required for the bipartite property to 
hold (see below). T-systems which are not bipartite can also be defined, and in that case, the matrix A is 
not symmetric. 

Example 2.3. The first example of such a T system is the one described in (|1.3| . In that case, we take 
the matrix A to be as follows: 

(2-4) 4'5=Vw'. 

where J a ,0 — C—2I is the incidence matrix of the Dynkin diagram associated with a simply-laced Lie algebra 
g. The coefficients q a are all set to be — 1. However, it is always possible to renormalize the variables so 
that q a = 1 in these cases [13] . and we use this approach here. 

In particular, if Q — A r , = <5 Q ,,3+i + 5 Q ,/3-i- This is the case we solve in this paper. 

We note that another example of generalized T-systems appeared in the context of preprojective algebras 
and the categorification program of The explicit connection was made in [4], Example 4.4. 
Finally, define the (possibly infinite) matrix P with entries 

(2-5) P^ = S a , (kj+i+^-i)- 

Then we can rewrite (|2.3|) as 

(2-6) T a<j . fc+ iT a j, fc _ 1 = \[ T^f q, \[ Y[{T a , , f;k ) A «''-" • 

a,j j' a' 

In the systems considered in 0], we allowed the matrix P to be a matrix with positive integer entries, such 
that it commutes with the matrix A, together with another condition on the sum of its entries (see Lemma 
12.51 below). Such a system is also a generalized bipartite T- system. 

2.3. Cluster algebra structure. We recall the formulation found in Appendix B of [3] of the cluster 
algebra associated with generalized (bipartite) T-systems. 

In the notations of Section [2j let S — (I r U I r ) x Z , and S = S U I' r . Each set I r , I r and I' r is just the 
set with r elements. For convenience, if a G I r , then by a we mean the ath element of I r , etc. 

We define the fundamental seed (x, B)q as follows. The variables Xo are 

x a ,j = T»j ; o, (a € I r ,j G Z); 

%, 3 = T at j. t i, (a ei r ,j e Z); 
(2.7) x a > = q a , a' € I' r ; 

(2.8) 

The elements of the set {x a j}U{xa,j} are the cluster variables and {a: Q '} are the coefficients. The exchange 
matrix of the fundamental seed is defined as follows: 

B a ,y,0,i = 0, (a, f3 G I r , j, I G Z), , = 0, (a,p G 7 r , j, I G Z), 

R _ — _ p3> 1 i Ali 1 — _ f>- 
(2.9) B a';0,j = = 

The last equation above denotes the entries of the extended /3-matrix, corresponding to the coefficients, 
which do not mutate. The matrices A, P are those of equation (|2.6[) for the generalized T-system. 
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Figure 2.1. A slice of the quiver graph of B, corresponding to constant a. The nodes 
in the strip are labeled by (j,k) of T aj; fc. The two subgraphs with even and odd j + k 
decouple in this slice, so we illustrate the only the connectivity of nodes of the same parity 
to node q a . The mutation /z reverses all arrows connected to q a . 



Example 2.4. In the case of the A r system (jl.lj) . we have the matrix A as in (|2.4|) . P as in (|2.5|) and 

q a = 1. In that case we do not need to include the coefficients q a , and the matrix B is equal to the matrix 
B. To recover the original T -system (jl.3[) . we take q a = — 1. 

It is clear that each of the mutations ^ a j and Ha.j exchanges one of the cluster variables in Xo via one 
of the T-system equation relations (|2 . 6[) . The mutation \i a .j acts on x as one of the T-system evolutions 
(|2.6p . where we specialize to k = 1: ^ a ,j{T a ,j;o) = T a j- 2 . Similarly, /i^jTaj-i — T QJ; _i is a T-system 
equation specialized to k = 0. 

Quite generally, if B a b = then jj, a o /i b = /i b o /i a . Since B a j-/3 t i = for all a,f3El r and j, Z € Z, when 
acting on the initial seed (x, B)o, the mutations \i a . m commute with each other for all a, m. Similarly the 
mutations Ha,m also commute among themselves. 

Therefore we can define the compound mutations 

which act on (x, B)q. More generally, Define (x, B)2k to be the seed with x a .j = T a j-2ki %a.j = 7a,i;2fc+i 
and B 2 k = B. Define (x, B) 2k+1 to be the seed with x a ,j = T a ,j ; 2k+2, Xa ,j = T a j. 2 k and B 2 k+i = —B. 
Then it is clear that fi(x 2 k) = ^-2k+i- Each mutation fi a j mutates the variable T a j- 2 k into the variable 
T Q j;2fe+2- Similarly, it is easy to check that 7t(x 2 fc) = x 2 fe-i, (i(3c 2k+ i) = x 2 fc and /Z(x 2 fe+i) = x 2fc+2 . 
The following statement is Lemma 4.6 of [3]: 

Lemma 2.5. Assume that the matrix A commutes with the matrix P , and that 
(2-10) E F ^ = 2 ^ 

k 

for any j. Then the cluster algebra X which includes the seed (x, B)o as in (|2.7[) . (|2.9|) includes all the 
solutions of the T-system (|2.6j) . All the T-system relations are exchange relations in this cluster algebra. 



To prove this Lemma, we need 
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Figure 2.2. The local action of the mutation \i on a section of the quiver graph. The 
compound mutation reverses all arrows connected to q a . 

Lemma 2.6. 

fi f(x, B) 2 kj = (x, B) 2 k+i, ~p ((x, B) 2 kj = (x, B) 2 k-\- 

Proof. In light of the preceding discussion, all that needs to be proved is that fi(B) = ~p(B) = —B. Let 
B' = /i(-B). Then, since -B Q j ; /3,fc = 0, we have 

• Ha,i{Bf3.j-^.k) = sigll(-B^ ,j;ct,i) [Bp ,j;ct,iB a ,i; 1 ,k\ + = 0; 

• /Jq,i(B^-. 7j( .) = -Bp jn k if (a,i) = (7, fc), and is otherwise unchanged, since if (a,i) ^ (7, k), 

Ma,i(%j ;7 ,fe) = ^,j; 7 ,fe + 8i « n ( B 3j;a,i)[- B 3,j ; a,i- B o.i;7.fc]+ = %j;7,fc- 

Similarly, (i a ,i(Bp t j?j t k) = -Bp J;ltk - 

• Recall the restriction that [P, A] = 0. Then 

• We have li a ,i(Bp<- a ,k) = —Bp>- 7i k, and otherwise, if (a,i) 7^ (7, fc) then /x Qi i 

^a,i(Bp'- ni k) = B/3' n ^k + ^ sign(B / 3' ;Qii )[i? / 3/ ;Qii i3 Q , ii;7j fc] = 5^ ;7 , 

so that /j,(Bp,. nt k) = -Bp,. lth . 

• Finally, using the restriction (|2.10p on the summation of elements of P, 

M(-B/3';7,fc) = + sign.(B/3>. a ,i)[Bf3>. atl B a ,i.^^k}+ = 8p n — ^ ^a,p ^ -Pq, 7 = ~ <5/?, 7 - 

In the quiver graph corresponding to -B, the last two statements are about how nodes x a j — T a ^ 2 k and 
Xa,j = T a j- 2 k+i are connected to node xpi — qp. If a ^ /3, they are not connected, and if a = (3, the 
connectivity is illustrated in Figure 12.31 and the mutations in Figure 12.31 
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We have shown that n{B) = —B. The proof that n{B) = —B is similar. 



□ 



Thus, we have shown that all the variables T a j.j~ appear in the cluster algebra, in fact, within a bipartite 
graph composed of the nodes reached from (x, B)q via combinations of the compound mutations /i and JL 
only. 

In this paper, we study the A r T-system solutions in terms of the fundamental seed cluster xo. The 
result will be an explicit interpretation of the solutions as partition functions of paths on a graph whose 
weights which are positive monomials in the variables xo. This will imply the positivity property [8] for 
the cluster variables T a .j-x- They can be expressed as Laurent polynomials with non-negative coefficients 
in terms of the initial data. 



3. Basic properties of the T-system 

From here on, we specialize the discussion to the T-system Note that the equation (jl.ip is a 

three-term recursion in the index k, and allows to determine all the {T a j t k+i}aei r ,jez m terms of the 
{{Taj^ki 2a,j',fe-i}aei'r,jeZ- We wish to first study the solution T a j fc to Equation fll.l[) in terms of the 
"fundamental" initial data x = (T a j , T at j t i) a ^x r ,jeZi that is, x . The techniques used in this section 
are a straightforward generalization of the methods used for the Q-system in [5] . We therefore present the 
proofs of the theorems in the Appendix, as they use standard techniques in the theory of determinants. 

3.1. Discrete Wronskians and conserved quantities. We can express the subset of variables {T a ; j t k : 
j, k G Z, a > 1} as polynomials of the variables in the set {Tij^ : j,k G Z}, cf [17]: 



Theorem 3.1. 

(3.1) 



T a ,j,k = det (Tij 

— a+b : k+a+b—a—l ) 

Ka.b<a 



a G I r , j, k G 



The proof of this theorem uses the standard Pliicker relations, and is similar to the case of the Q-system. 
We therefore present the details of the proof in the Appendix, Section IA.21 

If we consider a = r + 1 in Equation (|3.ip . since T r+ \j : k = 1, we have the polynomial relation among 
the variables {Tij-k}: 



(3.2) 



l,j,k—r 



l,j+l,/c+l-r 



Ti 



j—X,k+X—r 
1 ,j,k+2-r 



\,j-r+l,k-X 



T 



r+2,)fc 



T 



l,j — r,k 
l,j-r+l,k+l 



TlJM+r-2 Ti,j^\ t k+r-l 



= 1 



T\ j+r-l,fc-l Ti J+r _2,fc 

This is the "equation of motion" for the system. Since ipj^ is a discrete Wronskian determinant, it remains 
constant for solutions of a difference equation. The difference equation can be found by taking the difference 
of two Wronskians and arguing that a non-trivial linear combination of its columns must vanish. 

Theorem 3.2. We have the following linear recursion relations 

r+l 

(3.3) T i,j-b,k+b(-l) h Cr+i-b(j - A) = j, k G Z 

b=0 

where the coefficients c r +i_b(j — fc) depend only on the difference j — k, with co(m) = c. r+ i(m) = 1 for all 
m G Z, and: 



(3.4) 



r+l 
a=0 



T lj+a , fe+a (-l) a d r+1 _ a (j + k) = j, k G Z 



where the coefficients d r +i- a (j + k) depend only on the sum j + k, with do(m) — d r +i(m) — 1 for all 
m G Z. 



s 
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Such linear recursion relations can be obtained by noting that T r +2,j,k — and expanding the corre- 
sponding Wronskian determinant along the first row or column. The key fact to be proven is that the 
minors depend only on the difference j — k or the sum j + k. The proof is presented in the Appendix, 
Section [O] 

By analogy with the case of the Q-systems [5j[6], we may still call the variables Cb(k) and (4(fc) integrals 
of motion of the T-system, as they depend on one less variable than T. Moreover, they can be expressed 
entirely in terms of the fundamental initial data for the T-system, Xo. 

Example 3.3. In the A\ case, we have 

T\,j,k - ci(j - k)Tij-x t k + i + Tij_ 2 ,fc+2 = 
Ti,j,k — di(j + k)Tij + i^+i + Tij +2 ,fc+2 = 

with the integrals of motion 



ciU) = 
di(j) = 







1 






-3,1 








-2,0 




-2,0 






1 




h Tl - 


f3,l 








f2,0 




f2,0 



An explicit expression for the conserved quantities of Theorem 13.21 is as Wronskian determinants with a 
"defect" : 

Lemma 3.4. The conserved quantities c m (j) (m = 0, 1, r + 1, j S Z) of Equation (|3.3[) are 

(3.5) C m (j) = k^<* x {Tl,j+n+a-b,n+a+b-2) 

l<6<r+2, _ fa^r+2-T7i 

for any n E Z. 

Again the proof uses the standard techniques, and is found in Section IA. 41 of the Appendix. 

4. Conserved quantities and hard particles 

4.1. Recursion relations for conserved quantities. The conserved quantities l|3.5p satisfy linear re- 
cursion relations, which allow us to express them in terms of the initial data xo. We use recursion relations 
on the size r, so we first relax the boundary conditions T r +ij k ~ 1 f° r ai l i> fc G Z. 
Consider the T-system: 

(4.1) t at j t k+lta,j,k-l = t a .j + l.kt a ,j-l,k + t a +l,j,kt a -l,j,k: toj t k — 1) (j, k G Z, a G Z>o). 

Solutions of this system are expressible in terms of the initial data (ia.j.Cb ta,j,i)aeZ >0 ,jeZ- By definition, 
T a ,j,k — ta,j,k if we impose the boundary condition t r +i t j it = 1 for all j, k G Z. 

The proof of Theorem 13.11 does not involve the boundary condition T r +\ y t k = 1, so the determinant 
expression for t a j.k still holds: 

(4.2) ta,j,k = det (tij+a-^k+a+b-a-l) , a > 1. 

l<a,b<a 

Define the Wronskians of size N with a defect in position N — m: 

(4-3) CN,m,j,k = !^ a e < N (tl,j+a-b,k+a+b-N-l) , 

where cjv,m,j,fc = if m > N or m < 0. 

Note that cjv,o,j,Jfc = Tv,j,fc, and CN,N,j,k — Tv,i-i.fc+i by Theorem l3.H If we impose the second boundary 
condition of the T-system on the t's, then c r +i_ m j^k = c m (j — k + r). 
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Lemma 4.1. The Wronskians with a defect c a , m ,j,k satisfy the following recursion relations: 

(4*4) tot— — X,k— 1 C<x— X,m,j,k C<x— 2,m— l,j,k— 1 ~t~ ta — l,j,k Ca — l,m ,j— 1 

(4-5) ^a — lj — l.fc C-a,m,j,k ^a,j~l,fc+l Ca-l ; ni-l,j,fc-l ^a— l,m ,J— 1 

/or a > 2 and m, j,k > 1. 

Proof. The first equation (14. 4p follows from the Desnanot-Jacobi relation (|A.3p . with iV = a, ij. = 1, 
*2 = a, ji = a — ?7i, j*2 = a, for the matrix M with entries M ai h = Tij +a _b_i^-+a+fc-a-i, a,b = 1, 2, a. 

The second equation (|4.5p follows from the Pliicker relation (|A.2|) . with N = a, and the N x (AT + 2) 
matrix P with entries P Qi i = 5 a , a , and P a ,fc = T\ ,j+ a -fc.fc+a+6-a-i for 6 = 2, 3, a + 2 and a = 1, 2, a, 
and by further picking oi = 1, 03 = 2, 6i = a + 2 — m, and &2 = ct + 2. □ 

Theorem 4.2. TVie Wronskians with a defect c a _ m j,k defined in Equation ()4.3j) are uniquely determined 
by the following recursion relation, for a > 2: 

^a— l,k— 1 ta — l,j,k Ca.m,j.k—1 — ^a — — l,fe— 1 ^a.j~l,k Ca — 1.7n — l,j + l,k—l 
(4.6) "^ta j k — i t a — i j f» Ca — i^m.j — i^k — 1 ~t~ ^a.j,k~l ta,j — l,k Ca — 2,m—l,j.k~l 

and the boundary conditions co, m j,fc = 0m,o, /or aZZ j. k £ Z and ci tTn j.k = o~ m ,Qti,j,k + o~m,i ii,j-i,fe+i /or 
aZZ m, /, A: £ Z. 

Proof. Using Equation (|4.4p . the second line in (|4. 6|) is equal to tQ.^fc-i^a-i.j-i.fc-iCQ-i^m.j,*:- Canceling 
the overall factor i Q -ij-i.fc-i, we must prove that 

(4-7) t a — \,j,k Ca,m,j,k — 1 [tot,j— l.k Ca—l,m—l,j+l,k — l ~t~ ta,j,k—l Ca — l.m,j.k) 0* 

Multiplying the l.h.s. of (|4.7p by ia.j+i.fc and using (|4.ip we have 

ta,j+l,kta — l,j,kCa,7n,j,k — l (^a,j,fc+l^Q,J,^-l ^a+l,j/,fc^a— l,j,fc)Ca — l,m— 1 ,J+1,A>— 1 

ta,j-t-l,kta,j,k — lCa — l,m,j,k 
— ta — l.j,k(ta,j+l,kCa,m.j,k — l ~t~ ta+l,j.kCa — 1,771— l,j+l,fe— l) 

ta,j,k — 1 (^a,j,fc+lCa — l,m— l,j+l,fe— 1 ~l~ ^a,j + l,k^OL — l,m.j,k) 
ta,j.k—l(ta — l.j,kCa,m,j + l,k ta,j.k-\-lCa — 1,771— l,j+l,fe— 1 ^a,j + l,&Ca — l,7n,j,fc) 

where we have simplified the third line by use of (|4.4p . and finally used (|4.5p . Equation (|4.6p follows. 

Equation (|4.6p is a three-term linear recursion relation in the variable a and therefore has a unique 
solution c given the initial conditions at a = 0, 1. Moreover, these initial conditions are identical to those 
for (|4.3p . hence this solution coincides with the definition <|4 . 3 j) for all a,m,j, k. □ 

Let us define: 

(4.8) C a+hm {j,k) = Ca + 1 ™ +k - a '\ a>0, m,j,keZ. 

~ba>-\-l.j-\-k — a.k 

These satisfy C a+lfi (j, k) = 1 and C a+1>a+1 (j, k) = t a+1 j +k - a ^ lik+1 /t a+ltj+k _ atk . The conserved quan- 
tities of the A r T-system are obtained by imposing the boundary condition t r+1 j k = 1, in which case: 
c m(j) = Cr+i,m(j, k) for any j £ Z, and m = 0, 1, 2, r + 1, independently of k £ Z. 

Corollary 4.3. TVie quantities C a +i, m (j,k) of eg. (|4.8p are i/ie solutions of the following linear recursion 
relation, for a > 1 : 

(4.9) C a+ i, m (j, fc) = C a>m (j -2,k) + j/2a+i(j - a, /c) C a>m -i(j, fc) + y 2a (j - a - 1, fc) C Q _i, m _i(j - 2, fc) 
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2r-2 2r 2r+1 





3 5 2r-1 

Figure 4.1. The graph G r , with 2r + 1 vertices labeled i = 1, 2, 2r + 1. 



wzi/i coefficients: 



2/2a+lU,fc) = — (a > 1) 

ta+l ,j+k,k ta,j+kM+l 

2/2a(j,fc) = — (a > 1) 



(4.10) 2/i(j,fc) = , 

tl ,j+k+l,k 

subject to the initial conditions Co, m (j, k) — S m fi and Ci jm (j, k) = 5 m ,o + ^m,l — 1> 
Example 4.4. FFe have the following first few values of C a+ i t . m (j, k): 



a = 0: 


Ci, (i,fe) 


= i 








Ci,i(i,fe) 


= yiti - 






a = l: 


c 2 ,o(i,fc) 


= l 








c 2 ,iC?,*) 


= 2/1 (i - 


3,fc)+2/ 2 (i-2 ) fc)+7/ 3 0'-l,fc) 






C 2 ,2(i,fc) 


= yi(j - 


l,k)y 3 (j-l,k) 




a = 2: 


c 3 ,o(i,fe) 


= l 








c 3 ,i(i,fe) 


= 2/1 0' ~ 


5,k) + y 2 (j-4,k) + y 3 (j -3,fc) + 


Vi{j - 




c 3 ,2(i,fe) 


= 2/1 (j - 


3, fc)2/ 3 (j - 3, fc) + y 4 (j - 3, fc)yi(j 


-3,k) 






+2/5 (j 


-2,fc)(y 1 (7-3,fc)+y 2 Cj-2,fc) + 


2/3 (j - 




c 3 ,3(i,fe) 


= 2/1 (j ~ 


l,k)y 3 (j-l,k)y 5 (j-2,k) 





Remark 4.5. F/ie Corollary \4-3\ allows to interpret the conserved quantities of the A r T -system as follows. 
From the recursion relation (|4.9|) . we deduce that C r +i !m (j, fc) is a homogeneous polynomial of the weights 
2/1)2/2) ■••)2/2r+i, themselves ratios of products of some t a ^, c 's with c only taking the values k and k + 1. If 
we impose t r +ij.fc = 1, we see £/ia£, as explained above, C r +i jm (j, fc) = c m (j) is independent of k. We may 
therefore write C r +i !m (j, k) — C r +i,m(j) 0), tte latter involving only Ta^^'s with c = 0, 1. These give r 
conservation laws for m = 1, 2, r. For r = 1, me /lave /or instance 

Tij+k : k 1 Fi j+fc^fc+i 



Fi j'+fc-i^+i 2i J j+fc_i,fc+iTi )i j-+fc_2,fc Fi J +fc_2,fc 

Fl j,q 1 Fl,j-3,1 



C 2 ,i(j,0) 



Fi j^ijFi j-2,0 7ij-2,o 



4.2. Hard particle interpretation. In this paper, we introduce a slightly generalized model of hard 
particles on a graph. 

4.2.1. Definition of the model. Let G r be the graph of Figure |4~T1 with vertices labeled as shown. When 
r = 1, G\ is just the chain with 3 vertices, and when r = Go is a single vertex. 
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To each vertex labeled i in G r , we associate a height function h, where 

i+1 



h(i) 



, (»>1), h(l) = 0. 



A configuration of hard particles on G r is a subset S of /2r+i such that i,j £ I implies that vertices i 
and j are not connected by an edge. We can think of the elements of / as the vertices occupied by particles. 
The set of all hard particle configurations of cardinality m on G r is called C m . There is a natural ordering 
on the set I 2 r+i, and in the generalized hard particle model we define in this paper, the set S is considered 
to be an ordered set. 

In general, a hard particle model on G r associates weights to the occupied vertices which depend on 
the vertex label, and possibly also on the total number of occupied particles. The corresponding partition 
function is the sum over all possible hard-particle configurations of the products of the occupied vertex 
weights. 

For the purpose of this work, we define the partition function for m hard particles as 

m 

(4.11) Zg'(j, k) = J2 II fi«C7 - 2(r + t - m) - 1 + h{u),k) 

see m t=i 

with the weights yi as in (|4.10[) and S = {i\, i m }- 

4.2.2. Conserved quantities as hard particle partition functions. We have the following. 

Theorem 4.6. The partition function Z^ a (j,k) (|4.11|) for m-hard particles on G a coincides with the 
quantity C a+ i tm (j,k) of (|4.8|) . 

Proof. Hard particle partition functions on G r satisfy a recursion relation in r. Fix m and consider the 
configuration of particles on vertices (2r + l,2r). There are 3 possible pairs of occupation numbers for 
these two neighboring vertices, (0,0), (1,0) and (0, 1), respectively contributing to the partition function: 

• (0,0) contributes Z? y r x (j - 2,k). 

• (1, 0) contributes y 2r +i(j - r, k) Z m r Si (j, k). 

• (0, 1) contributes y 2r (j - r-1, k) Z m r S* (j - 2, fc). 
This implies that Z^ satisfies the recursion relation 

(4.12) Z%+* (j, fc) = Z%{j - 2, fc) + y 2r+1 (j - r, fc) Zg^j, fc) + y 2r (j - r - 1, fc) Z^U - 2, fc). 

But this is the same relation satisfied by C rjTn (j, fc), Equation (|4.9|) . with the same initial conditions, 
Zo r (j, fc) = 1 (for any r) and and Zf°(j, fc) = yi(j - 1, fc) = C 1A (j, fc). The theorem follows. □ 

Setting t r +i,j,k = 1, we have: 

Corollary 4.7. The conserved quantities c m (j) of the A r T-system are the partition functions for m-hard 
particles on G r , with the weights: 

/„• 7 \ T a+1 ]+k _ 1 k+1 T a j+1+k k 
y2a+iU) k ) = — ttt^ ^ — (1 < a < r) 

1 a+l,j+h,k J-a,j+k,k+l 

, ■ n T a+ l j + k k + l r a _l j+k+1 ,k n ^ ^ s 

V2 a {3, k) = = (1 < a < r) 

-La,j+k,k -ta,j'+fc+l,fc+l 



(4-13) Vl (j,k) 



T\ j+fe+i,fe 

where Toj t k = 2r+i,j,fc = 1 f or a ^ 3i k G Z. 

As the resulting hard-particle partition functions are independent of fc, we may set fc = in the expression 
for the weights. 
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h(v) 




FIGURE 4.2. A graphical interpretation of a hard-particle configuration on the graph Ge, with 
m — 5 particles at positions {1, 3, 6, 11, 13}. The label of the occupied vertex is indicated in a 
circle, and the time and height coordinates in rectangles. A distinct diagonal stripe corresponds 
to each particle. The leftmost stripe has x-intercept j — 2(r + 1) = j — 14, and the rightmost is 
at j — 2(r + 1 — m) = j — 4. The weight of this configuration is y 1 (j — 5, k)y^(j — 5, k)y&(j — 
6,fc)yn(j - 5,k)y rA (j - 6,k). 

4.3. A pictorial representation for the hard particle partition function. The hard particle con- 
figurations which give rise to the partition function of the form (|4.11|) can be represented graphically as in 
Figure S21 

• A particle at a spine vertex v (v G {1, 2, 4, 6, 2r — 2, 2r, 2r + 1}) is represented by a diamond on 
the two-dimensional lattice, its center at the height of the vertex, at the point {t, h(v)) for some 
t £ Z, and its vertices at the four neighboring lattice sites. 

• A particle a vertex v G {3, 5, 7, 2r — 1} is represented by the lower half of such a diamond. 

We call t the time coordinate, and h(v) the height. Each polygon is at (t, h(v)) contained in a diagonal 
stripe s, bordered by the lines y = x — (t + 1 — h(v)) and y = x — (t — 1 — h(v)). We denote s by its 
x-intercepts, s = {t — 1 — h(v), t + 1 — h(v)}. 

Given a configuration S G C m , with S = {i± < %2 < • ■ ■ < i m }, the polygon representing the particle 
i\ is drawn in the stripe s\ = {t — 2,i}; that of ii in the stripe immediately above and to the left, 
Si = {t — 4, t — 2}, and the fc-th polygon representing ik lies in stripe Sk = {t — 2k, t — 2k + 2}. The height 
of each polygon is determined by h(ij) and its time coordinate by its stripe: = h{ik) — 1 + t — 2{k — 1). 

If we choose t = j — 2(r + 1 — m), then Equation (I4.11|) can be written as 

m 

(444) Z^(j,k)= J2 n^®) 

see m e=i 

5. Path formulation and positivity 

We now give an expression for T a j^ as a function of the initial data Xo = (TJg^o, 2^9 ! j ) x)^eJ r ,jez- It can 
be interpreted as the partition function of weighted paths on a certain graph, with time-dependent weights. 
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1 



3 
-Q- 



4 



r-1 



r 



4' 



r-1' 



r+1 



r+2 



Figure 5.1. The graph G r , with 2r + 2 vertices. 



G 3 




1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 



FIGURE 5.2. The planar representation of a typical path in 7 \ & on the graph G3. 



That is, we generalize the notion of a weighted path, so that the weight of a step in the path depends on 
the time at which it is taken. 

As a corollary of the formulation in this section, we have the positivity Theorem I5.7I for the variables 
T a ,j,k as a function of the initial data. 

5.1. Definitions. Let G r be the graph in Figure 15.11 It has 2r + 2 vertices, which are ordered as 
(0, 1, 2, 2', 3, 3', ...r, r' , r + 1, r + 2). Its incidence matrix A is 

A m , m ' = A m <. m = 1, (2 < m < r); A m . m+1 = A m+ i tTn = 1, (0 < m < r + 1). 

The vertex labelled is called the origin of the graph. We call the vertices i the spine vertices of G r , and 
the edges which connect i — > 1 ± 1 spine edges. 

We consider the set l ?^ b t2 of paths p on the graph G r , starting at time t\ and vertex a, and ending 
at time £2 > t\ at vertex b. We take ti £ Z, and each step takes one time unit. The path p may be 
represented by the succession of visited vertices, p = {p(t)) t =t 1 ,t 1 +i,...,t 2 > with p{t\) = a and pfa) = b and 
Ap( s ) lP ( s+ i) = 1 for any s. 

Let Wij (t) be the weight of a step vertex i to vertex j at time t. We define the weight of a path p S J'° 1 ' & t 2 
to be 

t 2 -i 

(5.1) w(p) = Y[ w p{s)Ms+1) (s), p(t x ) = a, p(t 2 ) = b. 

S— 1\ 

The partition function for weighted paths in J ) tft 2 ^ s 
(5-2) Z?l%= J2 

For later use, we define Z^f t2 — if t\ > t 2 . 

Paths can be represented on the lattice 1? as in Figure I5T2"1 We associate a vertical coordinate h(i) — 
h(i') = i to each vertex of G r . The horizontal axis is the time. A step a — > b at time t on G r is a step 
(t,h(a))^(t + l,h(b)). 
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S p S p 




(a) (b) 

FIGURE 5.3. The involution if between pairs of m- hard particle configurations and paths in 
^^(r+i-m) j+2k- Case (a): the first stripe traversed by p is absorbed into S" , which has m + 1 
particles, and p' € 3 3 ^2(r+i-m)+2 j+2fc' Case 3 )'- The bottom stripe of S is absorbed into p' , 
now in 7° j '° 2(r+1 _ m) _ 2 j+2k , while 5" has only m- 1 particles. 



We claim (see Theorem l5.5[) that there exists a choice of weights w a ,b(s), as functions of xq = (T a jfi, T a ,j t i) a ei r .jez, 



such that Tij.k/Ti.j+kfl is equal to the partition function Z>j'_ k j +k . 

Dividing a path, which takes place from time t to time t' , into a first part from t to t', and a second 
part, from t' to i", we have 

h 

xeS T 



(5.3) z# = V Z a t ;*Z%„ t'e[t,t"] 



In particular, the matrix of one-step partition functions is called the transfer matrix 7(t), with 

entries 

(5-4) {7(t)) a , b = Z a t b t+l = «; 0|6 (t) A a<b . 

The transfer matrix is a decorated adjacency matrix. The recursion relation (|5.3[) implies 
(5-5) Z£* 2 = (7(ti)7(ti + 1) • ■ • T(ta - 1)) 0>6 . 

We use the following definition for weights w a ,b(s) of paths on G r : 

Wm,m'(s) = 1, ^m',m(s) = J/2m+l(s,0), (m £ {1, r}), 

(5.6) iw m>rn+ i(s) = 1, w m+ i !m (s) = y 2m {s,0), (m G {1, ...,r - 1}), 

«->o,i( s ) = !) ^i.o(s) = Vl(s,0), 
in terms of the weights yi(s, k) of Equation (|4.13[> . 

5.2. An involution on pairs of weights. We define an involution <p on the set C m x J , J !_ 2 ( r +i_ m ) j+2fe' 
consisting of hard-particle configurations on G r and paths on G r . 

Let (5,p) G C m x 3 , J -!_ 2 , r+1 _ m N J . +2fc , with m G {0, ...,r+ 1}, fc G Z + . We refer to the graphical 
representations of Figures EL"2l and IOI and we draw S and p on the same lattice (see Figure [53)1 . where 
S is represented between the diagonal lines y — x — (j — 2(r + l)) and y = x — (j — 2(r + 1 — m)), and p 
starts at (j — 2(r + 1 — m), 0), the x-intercept of the bottom stripe of S 1 . 
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The path p has an initial section po within the diagonal stripe {j — 2(r + 1 — m),j — 2(r — rn)}, 
consisting of u consecutive up steps and (i) a down step (s, u) — > (s + 1, u — 1) or (ii) two horizontal steps 
(s, u) — > (s + 1, it) — > (s + 2, u), where s = j — 2(r + 1 — m) + u. p \ po is then to the right of this initial 
stripe. 

Let a be a map from path steps of type (i) or (ii) on G r to the vertex set of G r . It is defined as follows: 

!2u - 2, 2 < u < r + 1; 
1, u = 1; 

2r + l, it = r + 2. 

cr((s,it) (s + 2,u)) = 2u-l. 

Remark 5.1. Graphically, steps of type (i) and (ii) in p can mapped precisely to the polygons representing 
particles on G r . A step of type (i) is the NE edge of a diamond (hence a particle on a spine vertex) and a 
step of type (ii) the upper edge of a half-diamond. The map a represents this correspondence. 

Denote by i the image of a step under the map a. We must now distinguish between two cases. 

• Case (a): If i < i x and S' := {i,ii,i 2 , -,U C C m+ i, define p' G 3 3 °f 2(r+1 _ TO)+2 j+2k to be the 
path with p'(j — 2(r + 1 — m) + 2 + x) = x for x = 0, 1, u (case (i)) or x = 0, 1, u — 1 (case 
(ii)), and p'(x) — p{x + 2) otherwise. 

• Case (b): If i > i\ or i m } ^ C TO +i, define S" = {i2,«3, ■■■,im} S C TO _i, the hard particle 
configuration with the right stripe removed. It is now drawn between the diagonal lines y = 
x — (j — 2(r + 1)) and y = x — (j — 2(r + 1 — (m — 1))). As for the path p', 

- If ii £ {3, 5, 2r - 1}, define p'(J - 2(r + 2 - m) + x) = x for x = 0, 1, h(ii), p'(j - 2(r + 
2 - m) + h(ii) + y) = h{ix) for y = 1, 2. 

— Otherwise, p'(j — 2(r + 2 — m) + x) — x for x = 0, 1, /i(ii) + 1, and p'(j — 2(r + 2 — m) + 
/i(ii) + 2) = h(ii). 

In both cases, p'(x) = p(x — 2) for the remaining times. 

Remark 5.2. Graphically the map can be visualized as follows. If the particle represented by po can be 
added to S while keeping the hard-particle condition, then we do this, while changing po so that it consists 
only of up steps, starting two steps to the right of the original starting point of p. Otherwise, perform the 
opposite operation, changing the first particle to a path segment. 

In view of the graphical description, the map ip is clearly an involution. Moreover it is weight-preserving: 
In Equation (|5.6j) . only the steps of type (i) or (ii) have a non-trivial weight. Moreover, w(tr(step)) = 
w(step) according to Equation (|4.13|) (setting k — 0). Therefore, w(S,p) = w(S)w(p) — w(S')w(p') = 
w{S',p'). We have 

Lemma 5.3. 

(5-7) E(-l) r+1 " m ^"(j,0)Z j l 2(r+1 _ m)j+2fc = 0, (j EZ,k> 0). 

m=0 

Proof. This is the partition function for pairs (S,p) € U^ Sm x -^j-2(r+i-m) j+2fe' wu ^ n an ex t ra factor 
(— l) r+1 ~ m wn ich ensures that the contributions of (S,p) and ip(S,p) cancel each other. □ 

We can also consider the sum in Equation (|5.7p in the case where k < 0. The sum is non-trivial in those 
cases only if k > — r — 1, since Z^f, = if t > t'. We extend the definition of ip: (p(S,p) = (S,p) if S = 
or if the path p has length zero and ii > 1. 

Lemma 5.4. 

(5-8) EVi) r+1 - m ^(j,o)^ 2(r+1 _ m)j _ 2i = i-iyz^u, o) 

m=0 
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G' 

where Z m r (j,k) is the partition function of m hard particles on G' r , the graph G r with vertex 1 removed 
(or the contribution to Z Gr in which vertex 1 is unoccupied). 

Proof. We apply the involution argument in the previous Lemma for k in the range — r — 1 < k < 0. Pairs 
(S, p) which are not invariant under tp cancel each other. We are left with the contribution of the invariant 
pairs. The latter always have p = and the vertex 1 unoccupied. □ 

Equations (|5.8|) are an expression for the initial conditions of the partition functions 2^1_ 2 ( r +i_ m ) j-2i 
with 1 < i < r — lin terms of hard-particle partition functions. 

5.3. The T-system solution Tij t k as a partition function of paths. Our main result in this section 
is the following. 

Theorem 5.5. 

(5.9) Ti j. k = T ltj+k ,o Z° fi k j+k . 

Proof. We will show that Sj.k — Ti t j+k,o ^j'-k j+k sa tisfies the linear recursion relation (|3.3p and coincides 
with Ti j k when k G {0, . . . , r} for any j G Z. Given that (|3.3[) has r + 1 terms, this implies Sj /. — Xij,fc 
for all other k. 
The sum 

r+l 

^ ( — l) m( V+l-m(j - k) Zf-k-2m,j+k 
m=0 

is equal to the sum in Equation (|5.7[) . since Z ( ^ l r {j 1 0) = c m (j). Therefore, it vanishes for all j G Z and 
k G Z+. This implies that Sj t k satisfies the same recursion relation (|3.3|) as Tij k. 
As for the initial conditions, we see from Equation (|5.8|) that S^fc satisfies 

i 

(5-10) ^ (-l) m Z^: m (j, 0) 5 , i _ m _ 2(r+ i_ i): „ i = ^ ''(j, 0) T hj _ 2 (r+i-i),o- 

m=0 

We will show that the variables Tij./c satisfy the same relations. Let 

(5.ii) wfio-) = £(-irzl* m (j,o) r y-=»(°+i-o- m , m > (0 <<<a + l) 

Using the recursion relations (|4. 12|) . we find 

W? a+1 (j) = Wf "(j - 2) + y 2Q+1 (j - a)W x G _l(j) + y 2a (j -a- l^TC? - 2 )- 

This is identical to the recursion relations (|4.12|) for Zf a {j,Q). Comparing the initial terms, we find 
that W^U) = 1 and Wf = Z G °(j,0) - %^Z o Go (j,0) = 0, so that W?°(j) = Z°°(j,0)\ yi(i _ l)=o . 
Moreover, 



wtrv) = i, 

- ^i Gl (i J o)-^^i^ G() (j,o) = 2 1 Gi (i,o)| !/l(j _3)=o, 

Jlj-2,0 

w 2 Gl (i) = Z 2 Gl (j,0)-%^iz i Gl (j,0) + %^Z Gl ( J ,0) = 0, 

where we have used the identity between Z and C, Example 14.41 the definitions (|4.13p , and T-system 
relations. In short, we have W Ga (j) = Zf a (j, 0)| B1= o, valid for all initial data a = 0, 1. This implies 

(5.12) W l G "(j) = Z^(j,0)\ yi=0 for all a. 
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Thus, W i " (j) is the partition function for i hard particles on G a with weight on vertex 1, or alternatively, 

vertex 1 unoccupied. Therefore, W i r (j) — Z i ' (j, 0). 

Therefore, Tij,fc and Sj,k satisfy the same recursion relation and have the same boundary conditions. 
The Theorem follows. □ 

The reasoning of Theorem 15. 51 can be carried through by considering paths with weights y a (j) replaced 
by the weights y a (j, k) of eq. (|4.10[) . We therefore have the following. 

Corollary 5.6. If we define the weights w a ^{s,p) as follows: 

iv m ,m'(s,p) = 1, w m > tm (s,p) = y 2m +i(s,p), (m G {1, r}), 
%,m+i(S)P) = 1) w m+1 , m (s,p) = y 2m (s,p), (ro G {1, ...,r - 1}), 
v>o,i(s,p) = 1, w 1<Q (s,p) = yi(s,p), 

with the yi{s,p) as in (|4. 1 3[) . then the following identity holds: 

(5.13) T lijtk =Tx !j+Ptk - p Z°'° pJ+p ({y a (s,k-p) ) j - p < s < j +p, 1 < a < 2r + 1}) . 

Proof. We may view this as a particular case of the translational invariance T a j & — ► T a j k+i of the T- 
system, namely that T a jf. is expressed as the same function of the initial data {Tgj f.—p, Tg j k- p +i} 8ei r j'ez 
as T aj:P is expressed in terms of {T^o, T^,]. }•£<=/,. j eZ . We deduce that Ti JiP = Tij+pfiZ^ !j+p ({y a (s, 0)}) 
has the same expression as Tij.fc = Ti.j^+fe-p in terms of y a (s,k —p), and the corollary follows, as the 
prefactor itself comes from the substitution T\ j+p t a — > T\^ +p t k—p- D 

5.4. General T-system solution 7 1 a ,j,fe : families of non- intersecting paths. We may interpret T^.j^ 
directly in terms of paths by use of the determinant expression of Theorem 13.11 for T a< j ^ in terms of the 
Ti t e t m- Indeed, given weighted paths on an acyclic graph T, say with partition function Z Ste for paths 
starting at vertex s and ending at vertex e, the Lindstrom-Gessel-Viennot formula gives an expression 
for the partition function of a non-intersecting paths on T (i.e. such that no to paths share a vertex) as 
Zs U :;s a ;e lt ...,e a = deti<i t j< a Z Siter We obtain: 

(5-14) = ^11/ I , / 2/, ., , ^°..., SQ;ei ,...,e Q 

where Z®f Sa . ei 6q stands for the partition function of families of a non- intersecting paths in the plane 
representation of Section [5.1i starting at the points s a = (j — k + 2a — a — 1,0), a — 1, 2, a and ending 
at the points &h — (j + k + a+ l — 2b, 0), b = 1,2, a. Alternatively, one may think of the partition 
function % Sl ,..., SQ ;ei,...,e Q as that of a "vicious" walkers (i.e. never meeting at a vertex) on G r , going from 
the root to the root, respectively starting at times j — k + 2a — a — 1 and ending at times j + k + a + l — 2a, 
a =1,2,..., a, each step corresponding to a unit of time. 

Interpreted in this way, the T a .j,k are manifestly positive Laurent polynomials of the initial data, via 
the weights yp(t) and the prefactor in (|5.14[) . We therefore have the: 

Theorem 5.7. The solution T a j t k of the T-system is expressed as a positive Laurent polynomial of the 
initial data xo = {r/3,j,0)^/3,j,i}/3e/r,jez for all a G and all j, k G Z. 

Proof. The statement is clear from the above discussion for k > a + 1 for which all the partition functions 
in the determinant (|5.14[) have the form Z t ' u with t < u, and therefore can be interpreted within the LGV 
framework. For k > however, from the structure of the T-system, it is clear that T a ^^ only depends on 
a finite part of the initial data {Tp ifi, Tpj^i, \f3 — a\ < k, \£ — j\ < k}. In particular, if < k < a+ 1, then 
as only the (3 > a + I — k > are involved, we may truncate the size of the T-system to some A r >, with 
r' = r — (a + 1 — k) < r. Upon renaming the initial data accordingly, we may interpret T a j t k as T a i j t k 
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in this new T-system, where a' = a — (a + 1 — k) = k — 1. For this a' > k — 1 the LGV formula applies, 
and positivity follows. Finally note that the expression of T a ,j,i-k hi terms of {TpjfijTpj^i} is the same 
as that of T a j t /. in terms of the reflected initial data {Tjg.^i, Tjg^o}, hence positivity follows for k < as 
well. □ 



6. Operator formulation and positivity in terms of mutated initial data 

Let A be the space of Laurent polynomials in the variables {Ta,j,fc}- We consider the invertible "shift 
operator" d acting on the infinite-dimensional vector space over A with basis {\t) : t £ Z}, withd|i) = \t— 1). 
It acts on the restricted dual space V*, with basis (t\ such that (t\t') — 5t,e, as (t\d = (t + 1|. We consider 
the algebra of formal Laurent series in d with coefficients in A acting on V. All operator relations which 
we derive below are considered in the weak sense, as identities between matrix elements. We also adopt 
the operator notation for diagonal operators in this basis, for example, w a ,b\t) = w a ,b(t)\t). 

6.1. An expression using operator continued fractions. Theorem 15.51 implies 

(6-1) ^i,3.fc = Tij +k>0 (7(j - k)7(j -k + l)---7(j + k- 1)) 

where the transfer matrix 7(s) is defined in Equation (|5.6| . 

We define the operator-valued transfer matrix T to be the matrix with entries (i|T a .b = 7 a ,b(t){t + 1|. 
We also define operator- valued weights Y Q , such that 

(6.2) (t\Y a =y a (t,0)(t + l\ 

where the y's are defined in (|4.13[) . 
Using these, we can write 

(6.3) T 1J>fe - T lj+fe , (j - k\ (T 2fc ) Q Q \j + k) = T 1J+M (j - k\ ((I - T) _1 ) Q Q \j + k), 

where (I — T) -1 := J2„>o(T) n . Therefore, the operator F = ^(7 — T) _1 ^j generates the variables Tij^. 

To compute F, we row-reduce the matrix I — T. The result can be written as a non-commutative 
continued fraction: 

(6.4) F = 

l-d(^l-d(l-dY 3 -d(---(l-dY 2r _ 1 -d{l-dY 2r+1 )- 1 Y 2r )- 1 ---y 1 Y^ Y^j . 

Alternatively, we can write F = Fo where the operators Fj are defined inductively: 
F r+2 - 0, 

F fe = (1-dYafc-i -d¥ k+1 Y 2 ky 1 , (fc = r + l,r,...,3,2), 
Fi = (l-d¥ 2 Y 2 )-\ F ^(l-dF 1 Y 1 )- 1 =F, 

where each term is understood formal power series in d. 

This expression is easily understood in terms of paths. Note that each time increment corresponds to 
an insertion of an operator d. An up step at height k followed by a down step contributes a weight dY 2 k , 
while a level step at height k contributes the weight efSi^fe-i- 

The operator generating function for paths above height k, F^, is obtained by shuffling the two following 
possibilities: (i) a level step pair k — > k — > k (ii) insertion of a path above height k + 1 between steps 
k — > k + 1 and k + 1 — > k (see Figure 16. ip . 



POSITIVITY OF THE T-SYSTEM CLUSTER ALGEBRA 



19 




d y 2k -i d d y 2k d 

FIGURE 6.1. The enumeration of paths on G r with time-dependent weights y a {t). The paths 
from and to height k which stay above height k (generated by F&) are related to those from 
height k + 1 by arranging any number of horizontal step-pairs k — > k — * k and up-down step-pairs 
k — > k + 1 — > k, in-between which we insert any path from and to height k + 1 staying above 
height k + 1. The operator weights are indicated on the bottom. 

6.2. Mutations and operator continued fractions. As with the Q-system of [5], we would like to 
have expressions for T a ,• & as functions of other possible initial data of cluster seeds in the cluster algebra. 
Cluster positivity means that they are positive Laurent polynomials in this data, and we can prove this 
by giving path generating functions on graphs with positive weights for them. The operator formulation 
introduced above was designed to allow us to do this in the case of special seeds of the form 

(6.5) x M = {T a ,j, ma +i\i = 0, l,a € I r ,j £ Z}, 

where M is a Motzkin path of length r: M = (mi, m r ) with |mj — TOj+i| < 1. 

This case is special, because it the non-commutative version of our construction in [5] for the Q-system. 
The only difference is that we must now use the operator-valued transfer matrix, instead of a scalar, to 
account for time-dependent weights. 

6.2.1. Compound mutations and restricted initial data. The cluster seeds in Equation (|6.5p are obtained 
from Xo by acting on it with a sequence of the compound mutations of the form 

(6-6) [i a = Y[ Va,j, Mq = II 

iez jez 

Note that the mutation matrix Bo has the property that B^ J a — if j ' ^ j' , hence fj, a ,j commutes with 
yitaj', so the compound mutations are well-defined. 

The mutations (|6.6|) act on initial data xo via the simultaneous use of all relations (jl.ip for all j G Z to 
transform T aj; fe_i — * T a j-k+i (forward mutation) or T a j-k+i — > T at j-k—i (backward mutation), the action 
being that of fi a when fc is odd and [la when k is even. Starting from the seed xo, and acting only with 
(|6.6[) generates a restricted set of cluster seeds. If, moreover, we require that each of the mutations be one 
of the T-system equations, we obtain only seeds of the form xm as in (|6.5|) . 

Remark 6.1. This is very similar to the situation of Reference [5], where seeds of type xm consist of 
variables {R a ,m a ; Ra,m a +i}- Here, we replace each variable R a , m with the infinite sequence (T , a j,m) i jez- 

6.2.2. Operator continued fraction rearrangements. In [5], we have shown that the generating function for 
Ri : n may be expressed in terms of any mutated seed xm via local rearrangements of the initial continued 
fraction in terms of the seed xo. Here, we give the non-commutative version of the starting point, which is 
operator version of the two rearrangement lemmas for fractions used in [S] . 
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The following are operator identities to be understood as identities between matrix elements of Laurent 
series in d. They are proved by a simple calculation. 

Lemma 6.2. Let A, IB be elements inA((d)). Then 

(6.7) 1 + d(l - Ad - M)- 1 A = (1 - d(l - M^A)' 1 
Lemma 6.3. Let A,B, C, U be elements inA((d)), with A + W invertible. Then 

(6.8) A+ (1 - d(l - HO^C)- 1 ! = (1 - U- d(l - dC'^B'^V 
where 

(6.9) A'=A + B, B' = CB(A + B)-\ C = d _1 CA(A + B)^ 1 ^ 

Remark 6.4. Lemma \6.2\ has a path interpretation. The r.h.s. of equation f|6 . T|) is the generating function 
for paths on the integer segment [0, 2], from vertex to 0, with operator-valued weights: 

w(Q -»• 1) = w(l 2) = d, iu(2^1)=B, tu(l->0) = A 

The l.h.s. of equation (|6.7p decomposes these paths into the trivial one (length 0, contribution I), and all 
the others, which start with a step — > 1 and end up with a step 1 — ► 0, with respective weights d and A 
(in this order). In-between, we have the generating function for "rerooted" paths, from vertex 1 to vertex 
1, which consist of arbitrary sequences of either steps 1 — > — > 1 (with weight Ad) or steps 1 — > 2 — ► 1 
(with weight dM). We call the rearrangement of this Lemma a "rerooting". 

6.2.3. General case: mutations as rearrangements. For each Motzkin path M, the solution of the T-system 
is can be expressed in terms of the initial data at xm as 

(6.10) Ti^ k = T hj+k _ miimi (j -k + m x |F M |j + k- m x ) 

for some operator continued fraction Fm- Our main claim is that this fraction is obtained from F (|6.4p via 
a succesion of applications of Lemmas 16.21 and 16.31 

For each Motzkin path we will define weights Yj(M), (i E fer+i), which are monomials in xm and d. 
The fraction Fm is a function of these. As in [5], we find that the effect of mutations on Fm >-> Fm' is the 
following: 

• If a = 1 , use Lemma 16.21 to write 

(6.11) F M = l + C ff M Y 1 (M). 

where (t|Yi(m) = Ti, t . mi +i/Ti : t+i. mi {t + 1|- Then Lemma T6.3I enables us to rewrite F M as Fm', 
a function of xm' • 

• If a > 1, apply Lemma \6. 31 to the part of Fm involving the weights Y j g(M) with (3 > 2a — 1. 
In both cases, the weights Y ( g(M) are transformed into weights Y / g(M / ). 

We will provide the precise construction and proof in Section 16.31 but for clarity, we refer the reader to 
Appendix [Bj where the example of is worked out completely. 

6.3. Paths on graphs with non-commutative weights. In this section, we define graphs with weights 
in A[d, d" 1 ]. The generating functions Fm are path partition functions on these graphs. The graphs are 
identical to those introduced in [5], and the weights contain exactly the same information contained in 
the two-dimensional representation of paths on these graphs introduced in [S] . The construction presented 
here is therefore a rephrasing of these paths in terms of operators. 

Remark 6.5. Although the two-dimensional representation of paths used here is identical to the one we 
used in [S], we did not, in the earlier paper, have use for the full information contained in this path 
representation. In particular, the horizontal coordinate ("time" in our language) had no interpretation in 
the context of Q- systems. Here, it corresponds to what is known as the spectral parameter in the T-system 
equations. 
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(a) (b) (c) 

FIGURE 6.2. The construction of the graph Tm (c) for some sample Motzkin path M (a) for 
Ag. All the pieces glued (b) are represented vertically. The skeleton edges in (c) are labelled 1 to 
19, the spine vertices to 13. 



Since one is used to reading lattice paths from left to right, we have chosen to act on the space V* 
instead of V. In that way, the order in which they act on the space is the same as the order the path is 
traversed. 

6.4. The target graphs T m . Let M be a Motzkin path. We decompose it into pieces which do not 
change direction: M = Mi UM2 U • • • , where Mj = (a,a + k,a + 2k, ■ ■ ■ , a + (li — l)k) with k — 0, 1 or — 1. 
The type of subpath is called k. All the graphs used below must be drawn vertically (see Fig |6.2[) . which 
makes unambiguous the notion of top and bottom edges. 
We construct a graph for each i as follows: 

• If k = then — G" , the graph Gi i of Figure 15.11 (represented vertically) , with its bottom and 
top edges removed. 

• If k = 1, then Tm; is a simple (vertical) chain with 2li vertices. 

• If k = — 1, then is the graph G'( (represented vertically) decorated with additional oriented 
"descending" edges b — > a with ij + l>6>a + l>l. 

We then glue the graphs: |rM i+ i is the graph obtained by identifying the top edge of Tm; with the 
bottom edge of Fm, +1 . Define T'^ = |Fm 2 I ' ' ' i an d Tm is together with one additional bottom 
and top edge and vertex. The graph Tm is rooted at its bottom vertex. 

Each graph thus constructed has a spine, namely a maximal vertical chain of vertices consisting of the 
unprimed vertices of the various pieces glued. We label the spine vertices consecutively starting from at 
the bottom (see Fig |6.2f c)). The vertices off the spine, which are attached only to a vertex i are labeled 
i'. We define the skeleton of Tm as the graph with all edges i — » j removed where i > j + 1. The edges of 
the skeleton, referred to as skeleton edges are labeled 1, 2, 2r + 1 from bottom to top (see Fig |6.2f c)). 
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A path traverses each edge of a graph in one direction or another, and in our formulation, we weight steps 
in each direction differently. Therefore, we now consider the non-oriented edges in Tm as doubly-oriented 
edges, each orientation corresponding to a different weight. 

Assign a weight Y aj b(M) to the edge a — > b. The weight of any edge away from the root, Yj^+i = 
Yjji — d. Skeleton edges Y Q (M) pointing towards the root are independent weights, which we define 
below. Weights on edges i — > i — k with k > 1 are defined as the following product involving only skeleton 
weights or their inverses: 

/i-k+i \ 
(6.12) Y M _ fc = J] Y a+ i i0 (Y B , a ;)- 1 (Y , iB )- 1 Y^fe+i,^. 

Va=i-1 / 

The ordered product is taken over edges from top to bottom, along a path from vertex i to % — k. Note 
that Y a ^ a ' — d. By inspection of (|6.12p we see that (£|Yj cx (t — k + 2|, hence Y^i-k "goes back in time" 
by k — 2 units. 

6.5. The positivity theorems for {T a j^}- We now write Tijj. as the partition function for paths on 
Tm ■ The values of the skeleton weights arc determined by considering the effect of a mutation on the seed 
data - They are determined by a recursion relation, which can be solved explicitly. 

6.5.1. Transfer matrices and mutations. As we illustrated in [5], any Motzkin path has a unique expression 
as a sequence of forward mutations, mp i— ► m'p = mp + Sfj. a where M = (mi, m r ) and M' = (m^, m' r ). 
We restrict the mutations to those which increase m a by +1 only in the following two cases: 

• Case (i): m Q _i = m a = m a+ i — 1, 

• Case (ii): m a -i = rn a — m Q+ i, 

(together with their boundary versions). This restricted set of mutations is sufficient to construct all 
Motzkin paths in the fundamental domain. 

The initial step in the induction is the Motzkin path Mo. The path interpretation on Tmo = G r was 
given in Section [5l The operator transfer matrix Tm = T and the operator generating function Fm = ^ 
are expressed entirely in terms of the d operator and the skeleton weights Y Q (Mo) = Y Q (|6.2[) . 

The inductive step is as follows. Given Tm and its operator weights, consider a forward mutation /i a or 
fj,a- M i ► M'. These have associated transfer matrices Tm and Tm' corresponding to the graphs Tm and 
Tm'- We compare the associated generating functions Fm = (I — Tm)o~o an( i ^M' = (I — TmOo.o usm g 
the row reduction process. Both are operator continued fractions, which differ locally due to the struction 
of the graphs. We find that the two operator continued fractions are equal to each other if and only if the 
weights of the graph Tm' are related to those of Tm as follows. 

Theorem 6.6. Let Y' = Y(M') and Y = Y(M), where M' = /i a (M) or /i S (M). Ifa^l, then, 

• Case (i): 

Y' 2q ,_ 1 = Y2o!_1+Y2q 

(6-13) Y' 2a = Y 2a+1 Y 2a (Y' 2a ^ 1 )- 1 

Y2a+1 = ^ 1 Y2q+iY2 q -i(Y 2q _ 1 ) 1 d 

• Case (ii): in addition to the previous, we have 

(6.14) Y' 2a+2 = d- 1 Y 2a+2 Y 2a - 1 (Y' 2a _ 1 )- 1 d, and Y^ = d~ x Y fi d, V/3 > 2a + 3 

If a = 1, we simply have to substitute Y% — > c? _1 Yid in the above formulas. 

Proof. The proof is by Gaussian elimination as in [5 . The case a = 1 is special, as it requires a rerooting 
of the generating function. The transformation of weights must be applied on FJ^ = (I — Tm)7i as in 
(|6.1ip . which induces the substitution Yi — > d~ x Yid. □ 
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Let 

(.D.lOJ A Q ,t jm — Tp; ; Ha,t,m — 

-^a,t+l,m -^a — l,t+l,m 

Corollary 6.7. TTie skeleton weights obeying the recursions of Theorem \6.(A subject to the initial condition 
([O]) are the operators Y /3 (M), acting as (rlY^M) = yp(M;t)(t + 1\, with: 

^a.t-}-m a — mi — l,m Q 



(6.16) lfta_i(M;t) 



A a -i,t 



+m a _i-mi,m a -i 



(6.17) y 2Q (M;t) = Ma+1 - t+m °~ mi ' m ° +1 x J °A n +i,7 + °^m;!!° +1 if^-^+i + i 

^a.,t+m a -mi,m a 1 otherwise 



(6.18) 



if m a = m a -i — 1 
otherwise 



Proof. By direct check of the recursion relations (|6.13H6.1^|) . □ 

Thus, we have two expressions for the generating function of T\j^, one in terms of the seed data xm and 
the other in terms of the seed data xm' • We call the transition between the two expressions a mutation: It 
acts on the graph Tm and on its weights. Alternatively, it acts on the operator continued fraction expresson 
for Fm as a rearrangement. 

6.5.2. Positivity ofTij^. We note that the weights (|6.16[) y a (M.;t) are positive Laurent monomials of the 
initial data at xm- We therefore have a positivity result: 

Theorem 6.8. Ti.j^+mt/Ti.j+k.rn! is the partition function for paths on the rooted graph Tm with the 
weights of Theorem \6. 7[ starting from the root at time j — k and ending at the root at time j + k. As such 
it is a positive Laurent polynomial of the mutated data at xm . 

6.5.3. General solution and strongly non-intersecting paths. We now turn to the expression of T a ^ t k hi 
terms of the mutated initial data xm- We will interpret the determinant formula Theorem 13 . 1 1 for T a j t f. a 
la Gessel-Viennot, in terms of the strongly non- intersecting paths on the graph Tm introduced in [5]. 

Let us briefly recall the two-dimensional TM-lattice paths used to represent paths on Tm, given in [5]. 
There are fundamentally three kinds of oriented edges in Tm.- the horizontal and vertical "skeleton" edges, 
and down-pointing long edges, with weights which depend on the skeleton weights. The steps taken along 
these edges on Tm are represented in 1? as follows (see Fig l6.3l for an illustration): 

• A skeleton step i — » i + e, e = ±1, at time t becomes the segment from (t, i) to (t + 1, i + e) 

• A skeleton step i — » i' or i' — > i at time t becomes the segment from (t, i) to [t + 1, i) 

• A long step j — > i, j > i + 1 at time t becomes the segment from (t,j) to (t + j — i — 2, i) 

Note that the increment of x-coordinate for each step coincides with the time shift we have associated with 
each step. Indeed, all steps advance by one unit of time, except the long ones, which go back in time by 
j — i — 2 < 0. The only difference with [S] is that we now attach time- dependent weights to the steps namely 
a weight y a ,b{t) for a step a — > b starting at time t (In the operator language, we have operator weights 
Y a: b that act as (t|Y aj b = y a ,b(t)(t + h\, where h is the time-shift of the corresponding step, h = +1 for all 
steps except the long ones, for which h = a — b — 2.). We conclude that this representation is perfectly 
adapted to our weighted paths, as the x-coordinate is nothing but the time-coordinate. 
Let us consider T a j,k as a function of the initial data at xm. Writing 

(6 19) a i3^ J r m \ d&t a +fr)k+ a +fr — a ~ l+rni 

l~[b=l 2l,J+fc+26-a,mi l<a : b<a Tij + k+2b-a-l,m 1 

Using Theorem 16.81 we may interpret Tij^ a+ b i k+a+b-a-i+m 1 /Ti,j+k+2b-a-i,m 1 as the partition function 
for TM-lattice paths from s a = (j — k + a + 1 — 2a, 0) to e& = (j + k + 2b — a — 1, 0). The determinant is 
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FIGURE 6.3. Two-dimensional lattice path representation on the graph Fm, M = (2,1,0) of 
the A3 case. We have indicated the "up" steps (i.e. away from the root) and the "down" steps 
(towards the root). 




FIGURE 6.4. A typical edge intersection of Fiyi-paths (a) and the result of the flipping operation 
on it (b). We have indicated the weights of the steps. The paths in (b) are said to be "too close" 
to each other. 



simply a signed sum of products of such path partition functions, corresponding in turn to the partition 
function for families of paths starting at {s a }a=i an d ending at {ef,}" =1 , with the usual weights times the 
signature of the permutation of endpoints induced by the configuration. 

In the standard Gessel-Viennot case, these signs produce the necessary cancellations to only leave us 
with the contribution of non-intersecting paths, namely families in which no two paths share a vertex. This 
is best proved by introducing a sign-reversing involution that pairs up and cancels all the unwanted terms 
in the expansion of the determinant. 
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In the case of TM-paths, the situation is more subtle, as paths may intersect without sharing a vertex. 
In [5], we have produced an involution, which allows to interpret an analogue of the determinant (|6.19j) 
as the partition function of strongly non-intersecting TM-lattice paths. This involution consists in flipping 
paths as follows. We consider the first intersection between two paths within a family. If the intersection is 
at a common vertex, we interchange the portions of paths before the intersection. If it is not at a common 
vertex, we flip the two paths as indicated in Fig l6.4[ by switching their beginnings until the crossing. 

We make then the following crucial observation: 

Lemma 6.9. In the generic flipping situation of Fig \6J\ the flipped pair of paths has the same (time- 
dependent) weight as the original one, up to the sign of the permutation of starting points, due to the 
following relation: 

Vh+v,v{ u + h- 2)y i+k+v s +v (u + i + k-2) = y h+Vii+v (u + h- 2)y i+k+v , v (u + i + k - 2) 

Proof. By direct application of the formula (|6.12|) for the long edge weights. □ 

The only invariant families under this involution are those where the paths do not lie "too close" to 
each-other, as otherwise they get cancelled by applying a flip. 

Therefore all the conclusions of [5] still hold in the present case, and we have: 

Theorem 6.10. T a> j, k+mi / Hb=i ^a,j+fc+2i>-a-i,roi * s the partition function for configurations of a strongly 
non-intersecting T-^-lattice paths, with the weights of Theorem \6. 7[ As such, T a j k+mi is a positive Laurent 
polynomial of the mutated data at xm . 

7. Conclusion 

The T-system equations are a special case of a non- commutative Q-system. That is, one can write a Q- 
system for non-commutative variables such that its matrix elements coincide with the T-system equations. 
One can think of the non-commutative Q-system as a "non-commutative" cluster algebra. Special cases 
of non-commutative cluster algebras have been considered in several contexts, for example the quantum 
cluster algebras of Berenstein and Zelevinsky [3J, or the more general recursion relation introduced by 
Kontsevich [14] (in rank 2), with similar Laurent properties, and which can be solved in some special 
"affine" cases using our methods 17. The main idea is that the path formulation seems to be particularly 
well adapted to the explicit solution of such problems, and makes Laurentness and positivity manifest. 
This will be discussed in a future publication. 

We should mention that there has been a great deal of interest in T-systems with more restrictive 
boundary conditions [5D1 Q~2] . We hope that the construction introduced in this paper will provide a simple 
way of treating such boundary conditions and their consequences. 

Appendix A. Discrete Wronskians 

A.l. Pliicker relations. Let P be an N x (7V + fc)-matrix. Let p 61 >•••>'><' | be the determinant of the matrix 
obtained by deleting the k columns b\,...,bk of P, times the signature of the permutation that reorders 
these column indices in increasing order. Then we have: 

fc 

^ 2) |p<»i,...,Ofc| |p6i,.-,6fc| _ y ' |p6 p ,a 3 ...,a»| |j>6i,...,6 p _i,oi,6 p +i,...,6fc| 

p=l 

for any choice of 2k columns oi, a k and b\, of P. In particular, when k — 2, we have 

^ 2) |pOl,02 | |p&l,&2 | _ |pi>l,£>2 | |pai,£»2 | _|_ |pf>2,Q2 | |pi>l,Ol | 

for any N x (N + 2) matrix P. 

Equation (|A.2[) implies the Desnanot-Jacobi relation. Let M be an N X N matrix, and let |M|, M\, 
\M^f^\ denote the determinants of M, the minor obtained by erasing row i and column j of M, and the 
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double minor obtained by erasing rows ii,l2 and columns Ji,j2 of M, respectively. Let 1 < i\ < 12 < N 
and 1 < ji < j2 < N, then 

(A.3) |M| |M£;£ I = |M£ 1 |M£ I -\M»\\M£\ 

It is easily obtained as a particular case of eq. (|A.2p . for a 2 = j\, b 2 = J2, P%ai — Pi,bi = ^i,i 2 j and 
M is the matrix P with columns ai and b\ erased. Indeed, one checks directly that: _p ai ' b i = 
\pa 2 M\ = \Ml$\, |P q i^| = \Mjf\, \P b ^\ = -|Mf I, \P a ^\ = \M$\, and |P 6 i.°=| = -|Mf|. 

A. 2. T-system as discrete Wronskians. Here present the proof of Theorem 13 . 1 1 which uses the relations 
in the previous subsection. 

Proof. Consider Equation (|A.3j) . Let N = a + 1, i\ = ji = 1, 12 = 32 = N and choose the matrix M with 
entries M a ^ = Tij +a -b.k+a+b-a-2 for a, b = 1, 2, a+1. We denote by W a +i j,k = \M\ the corresponding 
"discrete Wronskian" determinant. Substituting this definition into eq. (|A.3|) . we have 

(A.4) Wa+i,j,fcW a _i, i , J: = W a j,k-iW a j, k +i - W a ,j-i } kW a j + i,k 

valid for a £ I r , provided we set Woj t k = 1- Note that Wij t k — Ti.j,k by definition. Comparing eq. (|A.4p 
with the T-system (jl.ip . we deduce that the T's and W's obey the same recursion relations and share the 
same initial conditions at a — and 1. As the system is a three-term recursion in a this determines the 
solution uniquely and therefore we have W a j,k = T a ,j,k for all a S I r , j, k G Z. □ 

A.3. Linear recursion relations. Here, we present a proof of Theorem 13. 2\ that the variables T a ,j;k 
satisfy linear recursion relations, with constant coefficients which are the conserved quantities. 

Proof. We perform the discrete analog of differentiating the Wronskian, and compute tpj^i^+i — fj.k = 0. 
Denoting by ipj^ = |gi,g2, --- igr+il and ipj-1^+1 — |fi , , • - ■ ,f r +i| as the determinants of column 
vectors gi, fi, we note that gi+i = ii for i — 1, 2, r. We may therefore rewrite 

(A. 5) tpj- 1:k +i - <fij,k = = |fi,f 2 , • • -f r ,f r+ i - (-l) r gi| 

with (f h ) a = Ti J _i +a _b ! fc +a _. r _i +b and (gi) Q = Ti iJ -_i +Q)fc+a _ T ._i. Therefore, there must exist a non-trivial 
linear combination of the columns of the matrix which vanishes. We write it as 

r 

(A.6) ^(-l) b Cj .+i- 6 (j, k)f b + co(j, k) (f r+1 - (-l)' r gl ) = 0. 

b=l 

Recall that the entries of the vectors ff, depend on j, k in a very particular way, namely (f&(j, k)) a +i = 
(tb{j + 1, k + l)) a , and similarly for gi. In order for the above linear combination to be non-trivial, we 
must therefore have Q,(j, k) = Cb(j — 1, k — 1) = • • • = q,(j — fc, 0) for all j, k e Z, hence the coefficients Cb 
only depend on the difference j — k. Finally, we may normalize the coefficients in such a way that Co = 1 
identically, and the first part of the Theorem follows. 

The second part is treated analogously, by considering the difference of Wronskians <^j+i,fc+i — <fj t k = 
and reasoning on the rows of the corresponding matrices. □ 

A.4. Conserved quantities as Wronskian determinants with defect. Here, we give the proof of 
Lemma 13.41 expressing the conserved quantities of the T-system as Wronskian determinants with defects. 

Proof. Let 7 m (j, n) denote the right hand side of of eq. (|3.5|) . 

It is clear that that 7o(j, n) — T r+ i J+ „.„ +r = 1 and 7 r +i(j, n) = T r +ij+ n _i, n + r +i = 1 as consequences 
of Theorem 13.11 and of the A r boundary condition. 

Let p € Z and define the (r + 2) x (r + 2) matrix D to be the matrix with entries — Ti.j +p+ i-b,p+b-i, 
and D a fi = T j+ n +a-b.n+a+b-2 for a = 2, 3, r + 2 and 6 = 1,2,..., r + 2. The identity (|3 . 3[) may be recast 
into a vanishing non-trivial linear combination of the columns of D, with coefficients c r+ 2-b(j)(~ l) b _1 , 
b = 1, 2, r + 2, hence the determinant of D vanishes. 
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Expanding the determinant along the first row, we find that 

r+2 r+2 

(A.7) = det(D) = ^(-l) b+1 A,6|^| = ^(-l) b -V+2-&0>)Tu+p+i-&, P +&-i, 

6=1 6=1 

as the determinants j m (j,n) are the minors \P[ +2 ~ m \. 

Since the Wronskian determinant T r +% j t k = 1 is non-zero, there exist no other non-trivial linear recursion 
relation than Equation (|3.3|) with strictly fewer terms, hence the coefficients in Equation (IA.7j) must be 
proportional to those in Equation (|3.3|) . As co(j) — "fo(j,n) = 1, we deduce that "f m (ji n ) — c m{j) for all 
m = 0, 1, 2..., r + 1, and the Lemma follows. □ 

Appendix B. Example of A 2 : rearrangements, graphs and paths 

Here, we illustrate the program of Section f6. 2. 31 in the case r — 2. We first present the rearrangements 
of the operator continued fraction F, which make positivity of R\_ n manifest in all three cases. Next, we 
interpret these in terms of partition functions for operator-weighted paths on graphs, to illustrate Section 



B.0.1. Rearrangements. The fundamental domain for the action of mutations on the fundamental seed xo 
is coded by the following three Motzkin paths with 2 vertices: mo = (0,0), mi = /ii(mo) = (1,0) and 
m 2 = /i2(mo) = (0, 1). We give below the three operator continued fractions corresponding to these points. 

Seed xo: The continued fraction -Fb(y) reads for the fundamental seed corresponding to the Motzkim 
path mo is: 



(B.l) 



F (y)= f 1 - d(l - d(l - dY s - d(l - dY6) _1 Y 4 ) l Y 3 ) 



with operators i e 4, acting as (t|Yj = Ui{t) (t + 1|, and: 

,.\ r,\ ^2,t,l U \ ^1,4+1, 0^2,4-1,1 

yiw = 7z , 2/2 (*) = ~ — ~ , j/3 w = — ^ — ~ 

-tl,4+l,0 -tl,4,0 J 1,4+1,1 J 1,4, 1-i 2,4,0 

/,n ^1,4+1,0 T 2 t+i,o 
2/4 (t) = , y 5 (t) = — 

J- 2,4,0-' 2,4+1,1 J 2,t,l 

Seed x 2 = /i 2 (x ): Following Section T6.2.3[ we apply the Lemma [6.31 to F , with A = Y 3 , B = Y 4 , 
C = Y 5 and U = 0: This yields F (y) = F 2 (w), where 




(B.2) F 2 (w) = ^1 -d^l -d(l -d(l -d(l -dW 5 ) _1 W 4 ) X W 3 ) 

with operators Wj, i G J5, acting as (i|Wj = Wi(t)(t + 1|, with: 

/.\ ,4,1 /.-i -^2,4,1 Ti t+1, 0-*2 t, 2 

™lW = 7f > W 2 {t) = , W 3 (t) = 

Jl, 4+1,0 J- 1,4,0 1 1,4+1,1 1 1,4, 1^2, 4+1,1 

/.\ T Xt+ ii ^24+11 

«"4(t) = ~ — ~ . m(t) = — 

-t 2,4,1-^2,4+1,2 J 2,4,2 

To obtain this, we have written W x = Yi, W 2 = Y 2 , and W 3 = Y 3 + Y 4 , while W 4 = Y 5 Y 4 W3 1 and 
W5 = d _1 Y5Y3W 3 " 1 d, and used the T-system to simplify the expressions. 

Seed xi = /ii(xo): Following Section f6. 2. 3[ we first apply the rerooting Lemma r6.2[ with A = Yi and 
B = VY 2 , where 

V = (1 - dY 3 - d(l - dYs) - ^^ -1 
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m, 





Figure B.l. The fundamental domain for A 2 , coded by the Motzkin paths mo,mi,m2, 
and the corresponding target graphs for the path interpretation, with their edge labels. 
Mutations are indicated by arrows. 



This allows to rewrite F (y) = 1 + di^ (y)Yj., with F '(y) = (1 - Y x d - dUY 2 )~ 1 . We may now apply the 
rearrangement Lemma EH with A = d^Yjd, B = Y 2 , C = Y 3 + (1 - dY 5 )- 1 Y 4 , and U = 0: this yields 
Fo(y) = 1 + dFi(z)Yi, where: 



(B.3) 



Fi(z) 



1 -di 1 -d(l -dZ 3 -d{l - dZ 5 ) _1 Z 4 ) (z 2 + (1 - dZ 5 ) _1 Z 6 ) 



Zi 



where Zi, i £ I5 act as (i|Zj = Zi(t)(t + 1|, with: 



Zl (t) 



Zi, 



Zi,t-i,iTj,t+i,i 



22,i-l,l22,i+l,l 
-il,t,l-il,i+l,2-t2,t,0 



*s(*) 



T\.t+l,\T2.t-2.1 
Tl.t, 2 T 2 t-lfi 



2,t,0 



To 



and Zg is a long step weight, expressed in terms of the skeleton weights as: Zq = Z 4 (dZ 3 ) 1 Z 2 (a particular 
case of Eq. (|6.12[) ). Note that it is diagonal, namely: 



(t\z 6 = z 6 (t)(t\, where z 6 (t) 



«a(t-l) 



1 



Zi i,aZ2,i— 1,0 



The above weights follow from the identifications: Z x = d~ 1 Y 1 d + Y 2 , Z 2 = Y 3 Y 2 Z7\ Z 6 = d- 1 Y 4 Y 2 Z 1 ~ 1 , 
Z 3 = d- 1 Y 3 d~ 1 Y 1 dZ^ 1 , Z 4 = d- 1 Y 4 d- 1 Y 1 dZ 1 ; 1 and Z 5 = d _1 Y 5 d, and the use of the T-system to simplify 
the expressions. 

B.0.2. Paths on graphs with operator weights. Recall first that the continued fractions Fo(y), F 2 (w) are 
such that Ti jjt/Ti j+f^o = (J — k\Fi\j+k), i = 0, 2, while, due to the re-rooting, we have = 
(j-k + llAij + k-l). 
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The three above operator continued fractions may be interpreted in terms of path counting as follows. 
We have represented in Figure [Bj] the three rooted target graphs r m , attached to the three Motzkin paths 
m = mo, mi, rri2, together with their edge labelings. We have the following 

Theorem B.l. For i = 0, 1,2, the quantities (t\Fi\t') are the partition functions for paths on the graphs 
r mi , from and to the root, starting at time t and ending at time t' , and with operator weights defined as 
the product over the operator weights for each successive step of the path, in the same order. The weights 
are d per step away form the root, and respectively Yi, 1{ and Wi per step towards the root, along the edge 
labeled i. 
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Proof. The proof is a straightforward adaptation of the argument of Section 16.11 it uses operator trans- 
fer matrices T<, and amounts to performing the Gaussian elimination of / — T,;, in order to compute 
((I — T,;) _1 ) , where indexes the root vertex on T mi . We give explicit expressions below. □ 

We now list the transfer matrices for the three cases above. In all cases, we have Fi = ((I 



To 



Remark B.2. Note that, as opposed to the two other cases, the transfer matrix Ti is not made of diagonal 
operators times d, as is diagonal, hence goes back one step in time compared to the other operators Zj, 
i — 1,2, 5. This necessity for the longer descending steps to go back in time was already observed in 5 
in the two-dimensional representation of the T m -paths. 
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